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Topological superconductors have been theoretically predicted as a new class of time-reversal- 
invariant superconductors which are fully gapped in the bulk but have protected gapless surface 
Andreev bound states. In this work, we provide a simple criterion that directly identifies this 
topological phase in odd-parity superconductors. We next propose a two-orbital U — V pairing 
model for the newly discovered superconductor Cua;Bi2Se3. Due to its peculiar three-dimensional 
Dirac band structure, we find that an inter-orbital triplet pairing with odd-parity is favored in a 
significant part of the phase diagram, and therefore gives rise to a topological superconductor phase. 
Finally we propose sharp experimental tests of such a pairing symmetry. 

PACS numbers: 



The search of topological phases of matter with time- 
reversal symmetry has been an active field in condensed 
matter physics[T]. In the last few years, a new phase 
of matter called topological insulators [2j |3] has been 
predicted t4j and soon experimentally observed in a num- 
ber of materials [5J [^. More recently, a new class of 
time-reversal-invariant (TRI) superconductors has been 
predicted by a topological classification of Bogoliubov-de 
Gennes (BdG) Hamiltoniansjll [8]. As a close cousin of 
topological insulators, the so-called "topological super- 
conductor" is fully gapped in the bulk but has gapless 
surface Andreev bound states hosting Bogoliubov quasi- 
particles[71 O HO]- Now the challenge is to theoretically 
propose candidate materials for this new phase. 

In this work, we first provide a simple criterion which 
can be directly used to establish the topological super- 
conductor phase in centra symmetric materials with odd- 
parity pairing symmetry. This criterion applies to super- 
conductors with spin-orbit coupling. We next study the 
possibility of odd-parity pairing in the newly discovered 
superconductor Cua;Bi2Se3[TT], which has a 3D Dirac 
band structure due to strong SOC. We propose a phe- 
nomenological model for Cua;Bi2Se3 with short-range in- 
teractions. Thanks to the peculiar Dirac band structure, 
we find a specific odd-parity triplet pairing is favored in a 
wide parameter range, giving rise to a topological super- 
conductor. We propose an unusual flux quantization in 
Josephson interferometry as a sharp test of such a pairing 
symmetry. We also explicitly demonstrate the existence 
of gapless surface Andreev bound states in the resulting 
topological phase. 

We start by introducing Nambu notation = 

[Ck,aQ' c^k,a/3(**!y)/3a]' where a, (3 =t, i label electron's 
spin and a labels the orbital basis for cell-periodic Bloch 
wave- functions. The BCS mean-field Hamiltonian H = 
dk £,l^'H(k)£,]i uniquely defines a BdG Hamiltonian 

n(k) =. [Hoik) - /i]r, + A(k)r„ (1) 

where Hq describes the band structure of normal metal, fi 
is chemical potential, and A is pairing potential. For TRI 



superconductors, 07i(k)O ^ = 7i(— k) where O = isyK 
is time reversal operation. 

The BdG Hamiltonian H(k) of a fully gapped su- 
perconductor, which describes Bogoliubov quasi-particle 
spectrum, formally resembles the Bloch Hamiltonian 
of an insulator. An important difference, however, is 
that 7i(k) has particle-hole symmetry inherited from 
the doubling of degrees of freedom in Nambu space: 
S7i(k)S = — 7Y(— k) with S = SyTyK. Because of this 
extra symmetry, Schnyder, Ryu, Furusaki and Ludwig[7] 
and Kitaev[S] have shown that 3D TRI superconductors 
are mathematically classified by an integer invariant n 
instead of Z2 invariants for insulators |2l [3]. Despite this 
difference, since 7i(k) belongs to a subset of TRI Hamil- 
tonians, we observe that v = n mod 2 is nothing but 
its own Z2 invariant as explicitly defined in Ref . [T5 . It 
then follows that v — 1 implies a nonzero n and is suf- 
ficient (though not necessary) to establish a topological 
superconductor phase. 

A powerful "parity criterion" has been advanced by 
Fu and Kane to evaluate v efficiently for materials with 
inversion symmetry [4]. This motivates us to study topo- 
logical superconductors in centrosymmetric materials, 
for which the pairing symmetry can be either even or 
odd under inversion. It follows from the explicit for- 
mula for n[7] that even-parity ones cannot be topolog- 
ical superconductors. In this work we focus on odd- 
parity superconductors satisfying PHq{Is.)P = ffo(— k) 
and PA(k)P — — A(— k), where P is inversion opera- 
tor. We now provide a simple criterion for odd-parity 
topological superconductors: 

Criterion: a fully gapped TRI superconductor with 
odd-parity pairing is a topological superconductor, if its 
Fermi surface encloses an odd number of TRI momenta 
in the Brillouin zone. 

A special case of this criterion has been proved [T2] 
for certain triplet superconductors in which Ho{\s.) = 
i?o(— k) and A(k) — — A(— k), i.e., inversion simplifies 
to an identity operator P ~ I. Here we generalize the 
proof to all odd-parity superconductors, as needed later. 
Proof: Since P'H(k)P ^ 7i(— k), the parity criterion of 



2 



Ref.|l] does not apply directly. Instead, because Ar^^ an- 
ticommutes with t^, 7i(k) satisfy: 



P7i(k)P = W(-k), P = Pt, 



(2) 



Since the operator P defined here satisfies P^ = 1 and 
[P, Q] — 0, P can be used in place of P as an inver- 
sion operator for odd-parity superconductors. The cor- 
responding parity criterion with P reads 



(-i)'^-n^2™(r„). 



(3) 



Here Fq's (a = 1, ...8) are eight TRI momenta in 3D Bril- 
louin zone satisfying Ta — —Ta up to a reciprocal lattice 
vector. £,2m{^a) = ±1 is the P eigenvalue of the 2m-th 
negative energy band at Ta, which shares the same value 
S.2miTa) = ^2m-i-i(rQ) as its Kramcrs degenerate partner. 
The product over m in ([s]) includes all negative energy 
bands of Ti{k). The physical meaning of ^ becomes 
transparent in weak-coupling superconductors, for which 
the pairing potential is a small perturbation to Hq. As 
long as the bands £„(ra) of Ho{Ta) stay away from the 
Fermi energy (which is generically true) , pairing-induced 
mixing between electrons and holes in the eigenstates 
i'mi^a) of H{Ta) Can be safely neglected. So we have 
6m(rQ) = P2miTa) X T2m{Ta) with p and T being the 
eigenvalues of P and separately, ([s]) then factorizes 
into two products over p and r. Now the key observa- 
tion is that the set of all negative energy eigenstates of 
H corresponds to the set of all energy bands of Hq (both 
above and below /x), which form a complete basis of Hq. 
So we find YlmP^miTa) = Det[P] = ±1 independent of 
Ta, and thus Yla.mP'^rniTa) = 1 . (jsl then simphfies to 



(-ir-n^2™(r„)=n(-i) 



(4) 



Here N{Ta) is defined as the number of unoccupied bands 
at Ta in the normal state. (|4| now has a simple geomet- 
rical interpretation: = or 1 if the Fermi surface of 
Hq encloses an even or odd number of TRI momenta, 
respective lyt2^]. The latter case corresponds to a topo- 
logical superconductor. 

A well-known example of odd-parity pairing is super- 
fluid He-3[14j. In particular, the TRI and fully-gapped 
P-phase has been recently identified as a topological su- 
perfluid [7] [3 IS], in agreement with the above criterion. 
This identification explains the topological origin of its 
gapless surface Andreev bound states theoretically pre- 
dicted before[r5 . Odd-parity pairing in superconductors 
is less well established. A famous example is Sr2Ru04, as 
shown by phase-sensitive tests of pairing symmetry [T^. 
However, the observed signatures of spontaneous time re- 
versal symmetry breaking [T^ seem to prevent Sr2Ru04 
from being a TRI topological superconductor. 

In the search for odd-parity superconductors, we turn 
our attention to the newly discovered superconductor 



Cua;Bi2Se3 — a doped semiconductor with low electron 
density and Tc = 3.8ivr[TT]. A most recent angle-resolved 
photoemission spectroscopy experiment [18] found that 
the dispersion near center of the Brillouin zone F 
strikingly resembles a massive 3D Dirac fermion, being 
quadratic near the band bottom, and linear at higher en- 
ergy. Upon doping with Cu, the Fermi energy moves into 
the conduction band, about 0.25eV above the band bot- 
tom in the "relativistic" linear regime |TH]. To the best 
of our knowledge, this is the first discovery of supercon- 
ductivity in a 3D Dirac material, which motivates us to 
study its pairing symmetry. 

The Dirac band structure in the parent compound 
Bi2Se3 originates from strong inter-band SOC and can 
be understood from k ■ p theory[Tn]. Since the (lowest) 
conduction and (highest) valence band at F have oppo- 
site parity, general symmetry considerations show that 
the k ■ p Hamiltonian i?o(k) to first order in k takes the 
form of a four-component Dirac Hamiltonian |20j : 



i7o(k) = toFq -I- w(fcxFi + kyT2) + v^^k^T^, 



(5) 



where F^'s (i = 0, . . . , 3) are 4x4 Dirac Gamma matrices. 
The four components arise from electron's orbital (u) and 
spin (s). As shown by first-principle calculations [5} fT9]. 
the conduction and valence bands of Bi2Se3 mainly con- 
sist of two orbitals: the top and bottom Se p^-orbital in 
the five-layer unit cell, each mixed with its neighboring Bi 
P2-orbital (z is along c axis). The two orbitals transform 
into each other under inversion and we label them by 
az = ±1. The Gamma matrices in Hq are then expressed 
as follows: {Fo,Fi,F2,F3} = {a^,az ^ Sy, -a^ 'S> s^, ay}. 

To study superconductivity in Cu2:Bi2Se3, we con- 
sider the following phenomenological effective Hamilto- 
nian with short-range density-density interactions: 



dx 



U 



n + 2Vnin2 



, (6) 



where na{x) = X]q=| j, '^aai^)^aa{x) is electron density in 
orbital a. U and V are intra-orbital and inter-orbital in- 
teractions, respectively. All other local interaction terms, 
such as (cVajc)^ and (c'axsc)'^, are neglected [23]. iJeff 
is to be thought of as an effective low-energy Hamilto- 
nian, which includes the effects of both Coulomb and 
electron-phonon interactions. We will assume that at 
least one of them is positive, giving rise to pairing. Since 
U and V are difficult to estimate from a microscopic the- 
ory, we will treat them as phenomenological parameters. 
Naively, one would expect that the intra-orbital effective 
phonon-mediated attraction would be stronger than the 
inter-orbital one. However, since the same is true for the 
Coulomb repulsion, it is possible that the overall effective 
interactions satisfy, e.g., < V < U. 

To determine the pairing symmetry of the U—V model, 
we take advantage of two facts: a) near Tc, A forms 
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TABLE I: Pairing potential in mean-field BdG Hamiltonian 
of f/ — 1/ model, and their transformation rules. 



an irreducible representation of crystal point group; b) 
the mean-field pairing potential is local in x and thus 
/c-independent. The form of all such pairing potentials 
A are listed in Table II, where r5 = rorir2r3 = GzSz 
and Fj-fc = iF^T/j. Also shown are transformation rules of 
A's under the following symmetry operations of Bi2Se3: 
inversion P = — Fq = —cTx, threefold rotation around 
the c axis C3 = exp(iFi27r/3) = exp(— is^Tr/S), and 
mirror about yz plane M = —iTi^ = —is^. We 
find that Ai,..., A4 correspond to Aig, Aiu, and i?„ 
representations of point group D^d respectively. The 
three A representations are one-dimensional so that 
the corresponding phases are non-degenerate. Among 



them, c Ai{iSy)c oc AiCi|Ci| + A[ci^C2i -f (1 ^ 2) 
is spin-singlet pairing with mixed intra- and inter- 
orbital (orbital triplet) components, which is invariant 
under all crystal symmetries; C"^A2(«Sy)c oc {ci^C2i + 
Cij,C2|) is inter-orbital (orbital singlet) spin-triplet pair- 



ing; c^A3{iSy)c (X 
spin-singlet pairing. 



(ci|Ci| — C2|C2x) is intra-orbital 
The Eu representation is two- 



dimensional with c'^A4(iSj,)c cx aci|C2| -l-/3ci|C2x, where 
a and (3 are arbitrary coefficients, leading to a SU{2) 
degenerate manifold at Tc- Of these phases, the A2 pair- 
ing phase is odd-parity, TRI, and fully gapped, with a 
Bugoliubov spectrum given by 



E±,k — 



\ 



^2) 



(7) 



where = y + (^^ + ^y) + ''^z^z- Since the Fermi 
surface only encloses the F point, according to our earlier 
criterion A2 pairing gives rise to a topological supercon- 
ductor phase in the U — V model for Cua;Bi2Se3. 

We now solve the linearized gap equation for Tc of the 
various pairing channels to obtain the phase diagram. 
For purely inter-obital pairing A2 and A4, the gap equa- 
tion reads Vx2a{Tc) = 1. For purely intra-orbital pairing 
A3, it reads UxsiTc) — 1. Here Xi{T) is the finite tem- 
perature superconducting susceptibility in pairing chan- 
nel Ai. A straight-forward calculation shows that 



X2 = Xo / dkd{eu - /x)Tr[F5oPk]V(2i^o). 



(8) 



Here xo = Dq detanh (^) /e, where Dq is density of 
states at Fermi energy and W is high-energy cutoff. The 




FIG. 1: Phase diagram of the U — V model, showing the high- 
est Tc phase as a function of m//x and U/V (assuming V > 0). 
The arrow shows the experimental estimate for m//i, which is 
about I |18j . Two phases Ai and A2 appear, which are even 
and odd under parity respectively (see Table |T|. The insets 
shows schematically the structure of the Cooper pair wave- 
function in the A2 phase, consisting of two electrons localized 
on the top (1) and bottom (2) of the five-atom unit cell. 



projection operator Pk = X]a=i 2 l0A,k) (</'A,k| is defined 
by the two degenerate Bloch states at k. As we will see, 
the integral over the Fermi surface in ([s]), which takes into 
account the interplay between pairing potential and band 
structure effects in a multi- orbital system, will play a key 
role in favoring A 2 pairing. The other two susceptibilities 
X3 and X4 can be obtained simply by replacing F50 in 
^ with F30 and Fio respectively. Using Pk = 5(1 + 
J2iy=o'^k^'^) ™^ = im,vkx,vky,Vzkz)/s]i for Dirac 
Hamiltonian Hq, we obtain X2 = Xo(l — w^Z/x^), X3 = 
X4 = 2x2/3. The gap equation for the intra- and inter- 
orbital mixed pairing Ai is: 



det 




UxoCi 
VX0C2 



0, 



(9) 



where C„ = {m/ixY^ for n = 1,2. From ^ and (|9|, 
we now deduce the phase diagram. Since xs < Xo and 
X4 < X2; ^3 and A4 always have a lower T^ than their 
counterparts Ai and A2, respectively. Only the latter 
two phases appear in the phase diagram. By equating 
their Tc's, we obtain the phase boundary: 



U/V = 1 - 2m^/n^ 



(10) 



Fig.l shows the highest Tc phase as a function oi U /V 
and m/ fi, for positive (attractive) V. The A2 pairing 
phase dominates in a significant part of the phase dia- 
gram. Note that experimentally, it has been estimated 
that m//i « I [HI- When V <Q the Ai phase is stable 
for U > rr? I \J}\y\, whereas for smaller \J the system is 
non-superconducting. The fact that the phase boundary 
starts at the point U = V and m = is not accidental: 
at this point the Hamiltonian (l6| has an enlarged [/(I) 
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FIG. 2: Phase sensitive experiments to test A2 pairing, which 
is odd under both inversion (r — > — r) and reflection about the 
yz plane (x — > —x). A superconducting ring made of either an 
s— wave (a) or d— wave (b) superconductor contains a segment 
of a A2 superconductor. The flux through the ring is nh/Ae 
(a) or {n + |)/i/2e (b), where n is an integer. 



chiral symmetry: c —^ exp(i9TQQ)c. Under the unitary 
transformation exp(j7rr5o/4), the two pairing potentials 

{isy)c and (FV^i^{iSy)c transform into each other[25 . 

From now on, we focus on the topologically non-trivial 
A2 phase. To obtain the surface Andreev bound state 
spectrum, we solve the BdG Hamiltonian in a semi- 
infinite geometry z < 0: 

+ v{k,Ti + kyT2). (11) 



The continuum Hamiltonian (11 1 must be supplemented 
with a boundary condition at z = 0. For a Cua:Bi2Se3 
crystal naturally cleaved in between two five-layer unit 
cells, the wave function amplitude on the bottom layer 
corresponding to cr^ = —1 must vanish, so that a^i) — 
'>p\z=o is satisfied. By solving H at kx = ky = 0, we 
find that a Kramers pair of zero-energy surface Andreev 
bound states il;± exist for fi^ > — Aj, i.e., as long as 
the bulk gap remains finite. The wavefunctions of ^± are 
particularly simple for to = [21^ : 



i'±{z) = e """^(cos fcozlcTz = 1) 
(g)|s2 = ±l,Ty = =Fl) 



sin kf)z\az 



-1» 
(12) 



where k — A2/VZ and ko — fi/vz- Using k • p theory, we 
obtain the low-energy Hamiltonian describing the dis- 
persion of surface Andreev bound states at small k^ and 
ky-. Hsf = Vs{kxSy — kySx)- The velocity Vs is given by 

vs = (V^-|«riT,|v+)/(V+IV'+) ^ vAl/^P. 

Finally, we discuss the experimental consequences of 
the A2 state. The topologically protected surface state 
can be detected by scanning tunneling microscopy. In ad- 
dition, the oddness of this state under parity and mirror 
symmetries has consequences for phase-sensitive exper- 
iments. Consider a c-axis Josephson junction between 
a A2 superconductor and an s— wave superconductor. 
Since the A2 state is odd under reflection about the yz 
plane, whereas the s— wave gap A^ is even, the lead- 
ing order Josephson coupling between the two supercon- 
ductors, — Ji(A*A2 -I- c.c), vanishes (as well as higher 



odd-order terms). The second order Josephson coupling, 
— J2[(A*)^A2 -I- c.c], can be non-zero (as well as higher 
even-order terms). Therefore, the flux through a super- 
conducting ring shown in Fig. 2a is quantized in units 
of 2^i22j. Alternatively, in a Josephson junction between 
a A2 superconductor and a wave superconductor ori- 
ented as shown in Fig. 2b, the first order Josephson 
coupling is non-zero. The flux through the ring in Fig. 
2b takes the value ^{n+}^) (n is an integer). The same 
holds for an s— wave superconductor which does not have 
the mirror symmetry relative to the yz plane. The ob- 
servation of these anomalous flux quantization relations 
would be a unique signature of the topological A2 state. 

To conclude, we present a theory of odd-parity topolog- 
ical superconductors and propose the newly discovered 
superconductor CuxBi2Se3 as a potential candidate for 
this new phase of matter. We hope this work will bridge 
the study of topological phases and unconventional su- 
perconductivity, as well as stimulate the search for both 
in centrosymmetric materials with spin-orbit coupling. 
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